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Distinct core and halo stellar populations and the formation history 
of the bright Coma cluster early-type galaxy NGC 4889 
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ABSTRACT 

We study the stellar population far into the halo of one of the two brightest galaxies in the 
Coma cluster, NGC 4889, based on deep medium resolution spectroscopy with FOCAS at the 
Subaru 8.2m telescope. We fit single stellar population models to the measured line-strength 
(Lick) indices (H/3, Mgb, [MgFe]' and <Fe>). Combining with literature data, we construct 
radial profiles of metallicity, [a/Fe] element abundance ratio and age for NGC 4889, from 
the center out to ~ 60 kpc (~ 4i? e ). We find evidence for different chemical and star for- 
mation histories for stars inside and outside 1.2R e = 18 kpc radius. The inner regions are 
characterized by a steep [Z/H] gradient and high [a/Fe] at ~2.5 times solar value. In the halo, 
between 18 and 60 kpc, the [Z/H] is near-solar with a shallow gradient, while [a/Fe] shows a 
strong negative gradient, reaching solar values at 60 kpc. We interpret these data in terms of 
different formation histories for both components. The data for the inner galaxy are consistent 
with a rapid, quasi-monolithic, dissipative merger origin at early redshifts, followed by one 
or at most a few dry mergers. Those for the halo argue for later accretion of stars from old 
systems with more extended star formation histories. The half-light radius of the inner com- 
ponent alone is estimated as ~ 6 kpc, suggesting a significantly smaller size of this galaxy 
in the past. This may be the local stellar population signature of the size evolution found for 
early-type galaxies from high-redshift observations. 

Key words: galaxies: halos — galaxies: individual (NGC 4889) — galaxies: abundances — 
galaxies: elliptical and lenticular, cD — galaxies: formation 



1 INTRODUCTION 

Brightest cluster galaxies (BCGs) are the largest and most luminous 
galaxies located in galaxy clusters. The merger, star formation and 
chemical enrichment history of these galaxies are imprinted in their 
kinematics and chemical abundances. Spatially extended measure- 
ments for such quantities can therefore constrain their evolution 
and formatio n, believed to be closely related to the formation of the 
cluster (e.g., |Pubinskilll998h and the presence of the di ffuse intra- 
cluster light (e.g.. lNapolitano et alj2003l : lMurante et al.ll2007t) . 

In the cold dark matter-based scenario, early-type galaxies 
(ETGs) for m via merging of subclumps with various masses. Ac- 
cording to lKobavashil (200J), their merging histories can vary be- 
tween spherical infall of gas-rich subunits at high redshift, and a 
sequence of merging events at different epochs and masses. Galax- 
ies of the former assembly history are created via a process that is 
similar to the classical monolithic collapse. This dissipative pro- 
cess gives rise to high metallicity in the galaxy center, and to 
significant logarithmic metallicity ([Z/H]) gradients that can be 
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steeper than —0-5 and correlate strongly with galaxy mass (e.g., 
IChiosi & Carrardl2002l : lKobavashill2004| - |Pipino et ai]|2008h . Gas- 
rich binary mergers, on the other hand, produce only shallow gra- 
dient s (~ —0.1) with weak mass dependence l lBekki & Shioval 
1999). Merging events that take place after most of the stars are 
in place reduce, but do not erase completely the pre-existing gra- 
dients, because of the way in which stars of diffe rent metallic- 
ities are redistribute d in both energy and radius JWhitdll98d 
di Matt eo et al . 2009). In case of a later gas rich merger with a cen- 
tral starburst, a metallicity gradient ca n form again, but is confined 
to the central regions after a few Gyr (Hop kins et alj2 009). 

Stellar population studies in ETGs based on spectral in- 
dices show metallicities higher than solar in their nuclei, negative 
logarithmic metallicity ([Z/H]) gradients ranging from —0.16 to 
— . 3 0, and nearly constant [a/Fe] element abundances with radius 
(e.g.,|K obavashi & Arimoto 1999; Sa nchez-Blazquez et aiT 2006b; 
iReda et alj|2007l : lAnnibali et al.ll2007l) . For BCGs, several studies 
provide infor mation on the chemical abundances reaching up to ~ 
20 kp c radius ^Carter et al.ll999l : lBrough et alj2007l : lLoubser et al.1 
120091) . No significant gradients of [a/Fe] are observed, as in the 
case of normal early-type galaxies, while the range of [Z/H] gradi- 
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ents is somewhat wider, with values from —0.20 down to —0.58. 
Stellar ages are determined to be mostly old, > 8 Gyr. 

In this paper we analyse Lick absorption line indices from 
deep spectra for NGC 4 889, one of the two central BCGs in the 
Coma cluster (see, e.g. Gerha rd et al.l l2007t) . reaching unprece- 
dented radii of ~ 60 kpc (~ 4 R e ) in its outer halo (Section 2). 
We then derive age, [Z/H] and [a/Fe] radial distributions using sin- 
gle stellar population models (Section 3). Finally, we discuss the 
implication of these stellar populations properties for the forma- 
tion history of the core and halo of this galaxy (Section^. In what 
follows we adopt a distance of D = 92.7 Mpc (NASA/IPAC Ex- 
tragalactic Database) a nd effective radius R e — 33'.'88 = 15.23 kpc 
Jjdrgensenet all 19951) for NGC 4889. 



2 LINE-STRENGTH INDEX PROFILES FOR NGC 4889 

Combining new data for the outer halo of NGC 4889 with literature 
data we first construct radial Lick index profiles extending from the 

center to ~ 60 kpc radius (~ 4R e ); see FigureQ] 

The outer halo measurements are taken from ICoccato et all 

j2010L paper I). These data cover the radius range ~ 7 kpc to 60 
kpc and are based on 8 hr (run 1, PA=81°7, <7i ns t — 76 km s _1 ) 
and 5.5 hr (run 2, PA=— 8°3, crjnst = 96 km s _1 ) exposures, long- 
slit spectra obtained with FOCAS@ Subaru on Mauna Kea. Typical 
signal to noise ratios range from ~60 to ~10 for the innermost and 
outermost spectra, respectively. 

The data analysis is described in Paper I and includes 
wavelength- and flux-calibration, convolution to the spectral res- 
olution of the Lick system, correction for line-of-sight velocity 
broadening and correction for the offset to the Lick system. The 
sky background is evaluated from both offset sky exposures and 
dark regions in the scientific frames, using the former to identify 
regions of the slit in the latter that are free from stellar light contin- 
uum. Then we used the sky spectra from these regions (~5f6 from 
the galaxy center in run 1, and ~2'.9 in run 2); they have the same 
continuum as in the offset sky, and ensure the best removal of the 
variable sky emission lines. After sky subtraction, residual unre- 
solved emission lines are still present in the Mgl bandpass (run 2) 
and H/3 blue pseudo continuum (run 1). They are possibly caused 
by internal reflection in the spectrograph and are removed via a 
Gaussian fit (Paper I). Scattered light is measured in regions of the 
CCDs not illuminated by the slit. It amounts to ~ 4.5-5.0 ADUs 
(~20% of sky), also in the offset sky images, and is nearly uniform 
across the detectors in both runs; it is subtracted off the science 
frames together with the sky background. 

On the final spectra, H/3, Mgi , Mg 2 , Mgfc, Fe5270 and Fe5335 
line-strength indices are measured (see lWorthev et al]| 19941) , from 
whi ch we determine <Fe> = (Fe5270 + F e5335)/2 and [MgFe]' 
= y/Mgb (0.72 ■ Fe5270 + 0.28 • Fe5335). Their errors are deter- 
mined by Monte Carlo simulations which accounted also for the 
errors in the radial velocities. Special care is taken to quantify pos- 
sible systematic errors caused by sky subtraction. Residual contri- 
butions of up to ±2% of the subtracted sky spectrum are added to 
the galaxy spectrum. Larger residuals are excluded because they 
would produce detectable line features in the spectra (see Paper I). 
New index values are then determined for these new spectra. The 
systematic deviations of the new values are found to be no larger 
than the error bars, except for Mg2, which we therefore do not use 
to derive stellar population parameters. 

For the inner regions of NGC 4889, we use spatially re- 
solved Lick index measurements along the major axis from 
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Figure 1. Lick absorption line indices f or NGC 4 889 as proje cted onto the 
galaxy major axis. Diamonds: data from Mehlert et al. 12000); Cors ini et alj 
12008), averaged in radial bins. Here the second panel gives unprimed 
[MgFe], and H/3 values are shown without (blue) and with offset correction 
(black); see text. Filled circles: outer halo data from Paper I, also averaged 
in radial bins. Second panel gives primed [MgFe]'; see Table [T] Squares 
and crosses: measurements before averaging (dark and light green are from 
both sides of major axis spectrum, red and pink from both sides of mi- 
nor ax is sp ectrum). Stars: Cent ral apertur e measurements, from|j0rgenserj 
1 19991 red),lMoore et alj fcooj , magenta). [Sanchez-Blazquez et alj J2006al 
green). lLoubser et alj <2009l blue). 



iMehlert et al.l d2OO0t) . and along the minor axes from lCorsini et alj 
j2008l) . These authors provide [MgFe]=VMg6 < Fe > instead of 
[MgF e]'. We als o include cent r al ape rt ure measurements from 
J0rgense 3 J 19991); Moor e et al. I d2002l) : I S anch ez-Blazquez et al.l 
( 2006ah; lLoubser et al.l d2009h . For each index, all halo and inner 
data are combined in a single radial profile (see Fig.[TJ, by project- 
ing the minor axis data points onto the major axis using constant 
ellipticity e = 0.34 teender et al .11 19891) . 

Average profiles are constructed by grouping and averaging 
the data points in radial bins. For the halo data, we use 5 radial bins 
centered at 7, 12, 18, 32 and 59 kpc from the galaxy center. The in- 
ner literature data are similarly grouped in 6 bins, with average radii 
0.14, 0.43, 0.84, 1.7, 3.7, and 8.7 kpc. Errors for the average values 
are determined from the errors in the individual points and from the 
scatter in the bins. This scatter derives partially from observational 
errors, but also from systematic differences between major and mi- 
nor axis and between both sides of the same long-slit spectrum. To 
take these systematic effects into account, we compute the x 2 of 
the weighted mean with respect to the error bars of the m^n indi- 
vidual points in each bin. In case it is larger than mbm, the error of 
the weighted mean is increased accordingly. 

Within the final errors, the inner and outer datasets agree in 
the region of overlap (~7-15 kpc), exc ept in H/3 where the more 
accurate newer data give lower values. Trager et"aT] 1 2008) show 
that the H/3 linestrengths measured by Meh lert et alj d2OO0h for sev- 
eral Coma galaxies including the second cD galaxy NGC 4874 are 
stronger by an average 0.38A than all other literature values, possi- 
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bly due to a too short slit leading to incorrect sky subtraction. For 
NGC 4889 we also find that the H/3 aperture value determined by 
Mehl ert et al ] j2003h is larger by 0.32A than the four aperture mea- 
surements shown in Fig.[T] while the metal indices agree within the 
errors. We therefore construct a second, offset-corre cted set of in- 
ner av erage points for H/3 only, using the data from Mehle rtet al.l 
J2000L corrected downward by 0.32A ) and lCorsini et al.| J2OO8I . un- 
corrected). These averaged radial profiles are shown in Fig.[T]along 
with the individual measurements. 



3 SINGLE STELLAR POPULATION (SSP) 
PARAMETERS: [Z/H], [a/FE] AND AGE 

In all bins of the radial average profiles in Fig. Q] we determine 
luminosity- weighted metallicity [Z/H], [a/Fe] abundance ratio, and 
population age from fitting single stellar population models of 
IThomas et afl fefXB) to the line-strength indices H/3, Mgfo, <Fe> 
and [MgFe]' or [MgFe]. Their grid of models covers age values 
from 0.1 to 15 Gyr, [Z/H] from -2.25 to 0.67 and [a/Fe] from 
—0.3 to 0.5. We interpolated this model grid so that we have steps 
of ~ 0.4 Gyr in age, ~ 0.018 dex in [Z/H] and ~ 0.02 in [a/Fe]. 
For one of the central data points and a few error bars we had to 
extrapolate to metallicities up to 0.9. 

For each radial bin, we determine the best fit model from 
the measured H/3, Mgb [MgFe]' and <Fe> values by minimiz- 
ing x 2 - Errors for [Z/H], [a/Fe] and ages are computed by means 
of Monte Carlo simulations. The results are shown in Figure[2]and 
are listed in Table Q] Logarithmic gradients of the SSP parame- 
ters are determined from standard linear regression fits in log-log- 
diagrams. For given parameter P (denoting either [Z/H], [a/Fe], or 
log 10 [age/Gyr]), its logarithmic gradient AP is defined through 
P(R) - P(R e ) = AP ■ log 10 (R/R e ). The derived stellar popu- 
lation gradients are shown in Fig. [2] 

3.1 Results 

Over the entire radial range from the centre to 60 kpc, the metal- 
licity and [a/Fe] profiles in Fig. [2] can no longer be described 
as single power-laws. Rather, they show a "break" at ~ 18 kpc 
(~1.2 R e ): in the inner regions, the logarithmic [Z/H] gradient 
is negative and steep: A[Z/H] = -0.35 ± 0.02 (-0.49 ± 0.03) 
using the corrected (uncorrected) inner H/3 data, with metallicity 
that decreases from ~5 times solar in the centre to ~0.8 solar at 
18 kpc. The [a/Fe] abundance ratio in the inner l.2R e is nearly 
constant at a value of 0.37, i.e., ~2.5 times solar. By contrast, in 
the halo outside ~ 1.2 R e , the logarithmic [Z/H] gradient flattens, 
A[Z/H] = —0.1 ± 0.2, and the [a/Fe] profile shows a steep neg- 
ative gradient, A[a/Fe]= —0.68 ± 0.23, reaching solar values at 
~60 kpc. No peculiar features corres ponding to this "b reak" are 
seen in the surface brightness profile JThomas et alj|2007l) . 

The robustness of the "break" in the abundance profiles to pos- 
sible systematic errors caused by sky subtraction has been quanti- 
fied as follows. New SSP models are computed using Lick indices 
derived from spectra including up to ±2% residual sky contribution 
(Section[2]l. The corresponding range of SSP parameters is shown 
by the yellow regions in Fig. [2] and approximately corresponds to 
the detection limit for the "break" in the [Z/H] radial profile. Resid- 
uals of (unrealistic) ±8% of the subtracted sky spectrum would be 
required to remove also the "break" in [a/Fe] (green lines in Fig. 

The inferred SSP ages show an almost bimodal pattern. The 



halo of NGC 4889 at radii > 18 kpc is uniformly old, with SSP ages 
ranging between 10 and 13 Gyr. Younger ages are determined for 
the inner regions of NGC 4889: using the original H/3 values results 
in average SSP ages between 2 and 8 Gyr, with the lowest values 
inferred in the central kpc. With the offset correction to the average 
literature H/3 scale, average ages in the inner regions increase to 
8 Gyr, but some of the aperture data still imply quite young ages 
(~ 3 Gyr; Fig. 0. 

Young inferred SSP ages for BC G centers are not uncommon 
dFisher et alJl 1991 IPipino et alj|2009h . However, such young SSP 
ages, if real for NGC 4889 and not caused by the heterogeneous 
nature of the inner measurements, do not necessarily indicate that 
the main bulk of stars are young but can also be due to a contamina- 
tion of a small percentage of young stars. In fact, a component with 
age 1 Gyr and 10% mass fraction can bring the inferred SSP age 
of an otherwise old population down to ~3 Gyr (Trage retalj2000l 
Table 8), similar to the lowest central values in NGC 4889. Con- 
versely, the inferred SSP [Z/H] and [a/Fe] of a composite stellar 
popula tion follow approxim ately its V-band luminosity weighted 
values JSerra & Trageill2007l) , i.e., are only moderately sensitive to 
the [Z/H] of a minor young component (correspondingly less than 
the difference between the black and blue points in the top panel of 
Fig-H}- Thus the steep [Z/H] gradient and high [a/Fe] characterize 
the dominant stellar population in the inner half of NGC 4889. 



4 DISCUSSION 

From measurements at large radii combined with previous litera- 
ture data, we find evidence for distinct stellar populations in the 
inner half and outer halo of the nearby BCG NGC 4889. These 
point to a rapid formation process for the core and a more extended 
accretion phase building the halo, as we now discuss. 

Within ~ 1.2R e = 18 kpc, all [a/Fe] measurements give high 
values, implying a flat profile at ~ 2.5 times solar abundance ratio. 
This implie s rapid conversion of the gas into stars. Using a formula 
derived by Tho mas et alj (J2005) from simple chemical evolution 
models would predict that most of the stars within 18 kpc radius 
were formed in less than 0.1 Gyr. This is about one dynamical time 
at 18 kpc. However, given the current unc ertainties about SNI a pro- 
genitors and delay time distributions (e.g. lPritchet et al.ll2008l) . this 
timescale of 0.1 Gyr is not very secure. At the same time, the com- 
bined measurements support a steep inner [Z/H] gradient (at least 
A[Z/H] ~ —0.3 dex), while the published SSP age values inside 
18 kpc show considerable scatter, including both old (10-13 Gyr) 
and intermediate age (3-6 Gyr) values. These data are consistent 
with a superposition of a dominant old population with a minor, 
younger component (see previous section). 

The rapid formation and steep [Z/H] gradient of the stars in 
the inner half of NGC 4889 are reminiscent of predictions of a 
quasi-monolithic dissipative collapse model, in which stars form 
in a rapid burst progressing from the outside in , while the gas 
collapsing to the center is continuously enri ched jCarlberjl 19841 : 
lArimoto & Y oshii 19871: IThomas et alj| 19991) . In this process, log- 
arithmic gradients of ~ —0.5 can be reached. The chemical prop- 
erties of the inner NGC 4889 population are consistent with a sce- 
nario in which a dominant, old population formed rapidly at high 
redshift i n quasi-monolithic , dissipative merger collapse, as de- 
scr ibed by Kob avashil J2004I) in the context of hierarchical models. 
Subsequently, several of such units could have been involved in one 
or at most a few dry mergers, such that the original steep [Z/H] gra- 
dients are only partially erased - a major dry merger between two 
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Figure 2. Single stellar population parameters and their gra dients in NGC 
4889. Diamonds sho w values derived from the data of Mehlert et al. ( 2000); 
ICorsini et al.l J2008l) averaged in radial bins, without (blue) and with offset 
correction (black); see Fig.[T] Filled circles show halo values derived from 
the data of Paper I also averaged in radial bins; see Table [T] Stars show 
values derived from the aperture measurements given in Fig.[T] The dashed 
lines show the logarithmic gradients in [Z/H], [ct/Fe], and age fitted to these 
data (in two radial ranges for the first two), with numerical values given on 
the figure. The shaded yellow regions indicate the range of SSP parameters 
that would result if systematic residual effects in the sky subtraction at the 
level of ^ ±2% were present. The green lines correspond to (unrealistic) 
±8% residuals (see Section [3?71 and paper I). 



ellipticals w ith identical gradients already reduces the gradient by 
a factor 0.6 ddi Matteo et alj|2~009l) - and such that at most a small 
fraction of more iron-enriched stars is added to the stellar popula- 
tion. 

By contrast, the halo of NGC 4889 at radii larger than 
1.2R e = 18 kpc is characterized by near-solar metallicities with 
shallow gradient, a steep [a/Fe] gradient reaching solar values at 60 
kpc radius, and old ages (9-13 Gyr, possibly decreasing towards the 
outermost radii). The old ages and lower but still enhanced [a/Fe] 
indicate that the stars in the halo were formed at early times, but 
over longer time scales than the stars in the core, out of g as that was 
already iron-enriched by exp loding type la supernovae l lMatteuccil 
ll994l ; lBekki & Shiovall 199% . 

While galactic winds powered by supernova and/or AGN 
feedback from the collapsing central galaxy could explain early 
truncation of the star formation in the halo due to removal of the 
gas, they cannot explain the longer star formation timescale for the 
halo regions with respect to the inner galaxy. On the other hand, the 
properties of the stellar population in the halo of NGC4889 are con- 
sistent with the i dea that these stars were accreted fr om shredded 
satellite galaxies jAbadi et all 20061 : iNaab et alj|2009h . 

These satellites could have formed in lower density regions 
and sustained their star formation histories for longer times, un- 
til infall into denser environments and interactions truncated star 
formation. The more extended star formation would enable SN la 
to contribute iron to the gas and allow the build-up of near-solar 



metallicities and a-enhancement. The merger of these satellites 
with NGC 4889 could have happened long after the star forma- 
tion was truncated. The total luminosity in the NGC 4889 halo 
is Lhalo ~ 10 1 Lb,q, integrating the surface photometry for 
R — 1.2 — 4i? e . This light is equivalent to a few Milky-Way size 
galaxies. The merging of a few such satellites with a BCG is ex- 
pected in merger trees compute d for BCGs in massive clusters like 
Coma (e.g. lMurante et al]2007h . 

The bi-modal nature of the stellar populations in the core 
(R < 18 kpc) and halo (18 kpc < R < 60 kpc) of NGC 4889 
is also consistent with recent resu lts on the size evolution of early- 
type galaxies with red shift (e.g iDaddi et ail 120051 ; iTruiillo et al.l 
l2007l : lcimatti et alj|2008l) . such that ETGs at z ~ 1 (z ~ 2) have 
sizes a factor of 2 (3-5) smaller than ETGs with sim i lar mass today 
(Ivan der Wei et alj 12008 1; Ivan Dok kum et alj [20081 ; ISaracco et all 
2009). Recently, van Dokk um et alj d2010h find from stacked rest- 
frame R-band images that massive ETGs have nearly constant 
mass within 5 kpc with redshift, but increase their envelope mass 
by a factor ~ 4 since z = 2, with effective radius evolving as 
R e oc (1 + z)' 1 - 3 . Scaling the present R e of NGC 4889 with this 
relation would predict R e ~ 6.2 kpc at z = 1. An estimate for 
the size of the inner [a/Fe] -enhanced population in this galaxy is 
the radius enclosing half of the current luminosity within 18 kpc, 
obtain ed by truncating the surface brightness profile ( Thom aTet alj 
12007b at that radius; this is also ~ 6 kpc. The consistency of these 
numbers suggests that we may have found local stellar population 
signatures of the observed ETG size evolution. 

In general, BCG galaxies have a range of [Z/H] gradients and 
[a/Fe] (see Introduction), thought to arise from the relative in- 
fluence of the ea r ly collapse component an d subsequent mergers 
teobavasrdl200l ; |De Lucia & Blaizotll2007l) . In NGC 4889, a rel- 
atively large (~50%) fraction of the galaxy light appears to have 
been involved in the early collapse. Because of the steep [Z/H] gra- 
dient and high [a/Fe] this galaxy may be a particularly good case 
for distinguishing the different core and halo populations. It will be 
important to obtain similarly extended [a/Fe] and [Z/H] profiles for 
a larger sample of bright ETGs and to connect their analysis with 
the properties of the high-redshift ETG population. 
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Table 1. Line strength indices and stellar population parameters for NGC 4889, averaged in radial bins. 
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2.92±0.08 


9 4 1 -§ , 
»-t_2.4 


0.47°_;° 3 07 


U.J8_ .04 






3.67 


1.66±0.12 


5.36±0.11 


3.77±0.07 


2.65±0.08 


8-6 3 _l.o 


0.4 o_« 07 


n 440.08 






8.74 


1.67±0.27 


4.66±0.30 


3.69±0.15 
[MgFe]' 


2.91±0.14 




^0.23 


n 22 - 06 






7 


1.32±0.12 


4.68±0.17 


3.48±0.12 


2.57±0.11 


> 12.2 


u - lu -0.07 


n sn - 04 

u " :5u -0.07 


4.7 


8" (1) 


12 


1.33±0.13 


4.48±0.12 


3.33±0.07 


2.47±0.08 


> 12.2 


u - ui -0.05 




1.9 


7" (6) 


18 


1.56±0.12 


4.26±0.14 


3.07±0.09 


2.24±0.09 


> 12.2 


-n nq° 04 

u - UM -0.04 


n 4n - 08 


0.9 


13" (4) 


32 


1.71±0.16 


3.93±0.16 


3.03±0.1O 


2.48±0.12 


10.80l_f 95 


-0 07 10 
u - u '-0.06 


0.23°J 7 10 


0.4 


30" (5) 


59 


1.77±0.24 


3.16±0.32 


2.79±0.17 


2.52±0.16 


9.70i-f 95 


-0 17 012 
u - 1 '-0.09 


0.05°^^ 


0.2 


32" (4) 



Notes: Upper part of the table based on data from Me hlert et al 1 J2000h ;l Corsin i et al 1 fe008l) : lower part based on data from Paper I. Values for literature data 
refer to the case in which the Mehlert et al. H/3 values are offset-corrected by 0.32A, and provide [MgFe] rather than [MgFe]' (see Section[2]for details). Col. 1: 
Mean distance of the radial bins from the galaxy center, projected along the major axis. Cols. 2-5: average values for indices and errors in the bins. Cols. 6-8: 
Corresponding best fit age, metallicity and Q-enhancement with errors. Col. 9: ratio between counts from galaxy and background (both sky and scattered light), 
collected from all pieces of the slit that contribute to that bin. These counts are obtained in the wavelength range 5000 A - 5150 A, which is free from intense 
spectral lines, and therefore representative of the galaxy and background continua. Col. 10: total radial extent covered by the bin. The number of slit portions 
(data points in Fig.[T} used for this bin is given in parentheses. 



Carlberg R. G., 1984, ApJ, 286, 416 

Carter D., Bridges T. J., Hau G. K. T„ 1999, MNRAS, 307, 131 

Chiosi C., Carraro G., 2002, MNRAS, 335, 335 

Cimatti A., Cassata R, Pozzetti L., Kurk J., Mignoli M., et al. 

2008, A&A, 482, 21 

Coccato L., Arnaboldi M., Gerhard O., Freeman K. C, Ven- 
timiglia G., Yasuda N., 2010, ArXiv:1006.1894 (paper I) 

Corsini E. M., Wegner G., Saglia R. R, Thomas J., Bender R., 
Thomas D., 2008, ApJS, 175, 462 

Daddi E., Renzini A., Pirzkal N., Cimatti A., Malhotra S., et al. 
2005, ApJ, 626, 680 

De Lucia G., Blaizot J., 2007, MNRAS, 375, 2 

di Matteo P., Pipino A., Lehnert M. D., Combes E, Semelin B., 

2009, A&A, 499, 427 
Dubinski J., 1998, ApJ, 502, 141 

Fisher D., Franx M., Illingworth G., 1995, ApJ, 448, 119 
Gerhard O., Arnaboldi M., Freeman K. C, Okamura S., 

Kashikawa N., Yasuda N„ 2007, A&A, 468, 815 
Hopkins P. F., Cox T. J., Dutta S. N., Hernquist L., Kormendy J., 

Lauer T. R„ 2009, ApJS, 181, 135 
J0rgensen I., 1999, MNRAS, 306, 607 

J0rgensen I., Franx M., Kjaergaard P., 1995, MNRAS, 273, 1097 

Kobayashi C, 2004, MNRAS, 347, 740 

Kobayashi C, Arimoto N., 1999, ApJ, 527, 573 

Loubser S. I., Sanchez-Blazquez P., Sansom A. E., Soechting 

I. K„ 2009, MNRAS, 398, 133 
Matteucci F, 1994, A&A, 288, 57 

Mehlert D., Saglia R. P., Bender R., Wegner G., 2000, A&AS, 
141, 449 

Mehlert D., Thomas D., Saglia R. P., Bender R., Wegner G., 2003, 
A&A, 407, 423 

Moore S. A. W., Lucey J. R., Kuntschner H., Colless M., 2002, 

MNRAS, 336, 382 
Murante G., Giovalli M., Gerhard O., Arnaboldi M., Borgani S., 

Dolag K., 2007, MNRAS, 377, 2 



Naab T, Johansson P. H., Ostriker J. P., 2009, ApJ, 699, L178 
Napolitano N. R., Pannella M., Arnaboldi M., Gerhard O., 

Aguerri J. A. L., et al. 2003, ApJ, 594, 172 
Pipino A., D'Ercole A., Matteucci F, 2008, A&A, 484, 679 
Pipino A., Kaviraj S., Bildfell C, Babul A., Hoekstra H., Silk J., 

2009, MNRAS, 395, 462 
Pritchet C. J., Howell D. A., Sullivan M., 2008, ApJ, 683, L25 
Reda F. M., Proctor R. N., Forbes D. A., Hau G. K. T, Larsen 

S. S., 2007, MNRAS, 377, 1772 
Sanchez-Blazquez P., Gorgas J., Cardiel N., Gonzalez J. J., 2006a, 

A&A, 457, 787 

Sanchez-Blazquez P., Gorgas J., Cardiel N., Gonzalez J. J., 2006b, 
A&A, 457, 809 

Saracco P., Longhetti M., Andreon S., 2009, MNRAS, 392, 718 
Serra P., Trager S. C, 2007, MNRAS, 374, 769 
Thomas D., Greggio L., Bender R., 1999, MNRAS, 302, 537 
Thomas D., Maraston C, Bender R., 2003, MNRAS, 339, 897 
Thomas D., Maraston C, Bender R., Mendes de Oliveira C, 2005, 
ApJ, 621, 673 

Thomas J., Saglia R. P., Bender R., Thomas D., Gebhardt K., et 

al. 2007, MNRAS, 382, 657 
Trager S. C, Faber S. M., Dressier A., 2008, MNRAS, 386, 715 
Trager S. C, Faber S. M., Worthey G., Gonzalez J. J., 2000, AJ, 

120, 165 

Trujillo I., Conselice C. J., Bundy K., Cooper M. C, Eisenhardt 

P., Ellis R. S., 2007, MNRAS, 382, 109 
van der Wei A., Holden B. P., Zirm A. W., Franx M., Rettura A., 

Illingworth G. D., Ford H. C, 2008, ApJ, 688, 48 
van Dokkum P. G., Franx M., Kriek M., Holden B., Illingworth 

G. D., et al. 2008, ApJ, 677, L5 
van Dokkum P. G., Whitaker K. E., Brammer G., Franx M., Kriek 

M., et al. 2010, ApJ, 709, 1018 
White S. D. M., 1980, MNRAS, 191, IP 

Worthey G., Faber S. M., Gonzalez J. J., Burstein D., 1994, ApJS, 
94, 687 



